Composite Dirac fermions in graphene 
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Generalizing the notion of composite fermions to the case of " pseudo-relativistic" Quantum Hall 
phenomena in graphene, we discuss a possible emergence of compressible states at the filling factors 
v = ±1/2, ±3/2. This analysis is further extended to the nearby incompressible states viewed as 
Integer Quantum Hall Effect of composite Dirac fermions, as well as those that might occur at 
v = 0, ±1 as a result of (pseudo)spin-singlet pairing between the latter. 



The traditional interest in Quantum Hall Effect has 
been rekindled by the recent experiments on mono- and 
double-layers of graphene where the interplay between 
unscreened Coulomb interactions and pseudo-relativistic 
kinematics of the Dirac quasiparticles has long been ex- 
pected to harbor a host of novel phenomenal. 

In graphene mono-layers, the Integer Quantum Hall 
Effect (IQHE) plateaus were found at the integer values 
<?xy = v = (An + 2), n = 0,±1,±2, . . ■£ (hereafter, we 
put h = e = c = 1, and measure all the conductivities 
in units of e 2 /h = l/2n). Elaborating on the earlier in- 
sight of Refsi, this observation was readily explained 3 
by treating the low-energy excitations in graphene as 
(pseudo)relativistic Dirac fermions with linear dispersion 
and speed vf ~ 10 6 to/s. These quasiparticles carry a 
physical spin s = 1/2 and possess an additional orbital 
("pseudo-spin" or "valley") quantum number (hereafter 
referred to as R and L) corresponding to the double de- 
generacy of the electronic Bloch states in graphene. 

The resulting SU(A) symmetry of the non- interacting 
Hamiltonian survives the long-range Coulomb interac- 
tions, although it gets broken in the presence of the Zee- 
man and various additional short-range (Hubbard-likc) 
interaction terms. A number of implications of this sym- 
metry which generalizes independent rotations in the spin 
and valley subspaces have been explored in several recent 
papers^. 

For one, by drawing a parallel with the previous stud- 
ies of spin-unpolarized double-layer Quantum Hall sys- 
tems, it was argued that, apart from the Dirac kinemat- 
ics, the situation in graphene is similar to that occurring 
in the double-layer systems in the limit of vanishing inter- 
layer tunneling. Thus, a graphene analog of the Quantum 
Hall Ferromagnet was predicted to occur, which type of 
strongly correlated states would manifest itself as addi- 
tional (interaction-induced) plateaus at all the integer 
filling factors^. 

In a recent experiment, additional plateaus were indeed 
observed at a xy = 0, ±1, . . thus suggesting a complete 
lifting of the spin and valley degeneracies at the lowest 
(n = 0) relativistic Landau level (LL). As an alterna- 
tive interpretation, it was pointed out® that the behavior 
reported in Ref£ can also be explained by invoking the 
"magnetic catalysis" scenario of Ref. 7 where, contrary to 
the predictions of Refsi*, the valley degeneracy gets lifted 
only at the th LL, in agreement with the data of Ref£. 

Still awaiting its observation, however, is a graphene 



counterpart of the Fractional Quantum Hall Effect 
(FQHE). By analogy with FQHE in the conven- 
tional ("non-relativistic") two-dimensional electron gas 
(2DEG) with parabolic quasiparticle dispersion, one 
might expect that its graphene analog can also be stud- 
ied by adapting the idea of statistical flux attachement 
to the case of the Dirac fermions. 

In the present Letter, we discuss such a procedure, 
thereby setting the stage for a systematic analysis of the 
(potentially, much richer than in the case of the conven- 
tional 2DEG) realm of FQHE phenomena in graphene. 

The flux attachment recipe would usually be applied 
to a half- (or, more generally, l/2p-, where p is an inte- 
ger) filled uppermost LL, while the rest of the system 
would be treated as an inert incompressible background. 
In a generic SU (TV)-invariant system, the Chern-Simons 
Lagrangian implementing a transformation from the orig- 
inal electrons to A-component composite Dirac fermions 
(CDFs) takes the form 
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Here go = 2i:e 2 / €$vf ~ 3 is the bare Coulomb cou- 
pling, the vector potential A 1 — (0,—By/2,Bx/2) rep- 
resents an external magnetic field, and the matrices 
7 Q = (l,d x ,(—l) a d y ) act in the space of spinors *S> a 
composed of the values of the CDF wave functions on 
the two sublattices of the bi-partite lattice of graphene. 

In a multi-component system, the statistical flux pro- 
vided by the Chern-Simons fields a l a can be attached in 
a number of different ways, the choice between which 
should ultimately be determined by the nature of the 
ground state in question. Accordingly, there exist dif- 
ferent choices of the integer- valued matrix K , the only 
condition imposed upon which is that the transformed 
CDFs retain their (mutual) fermionic statistics. This re- 
quirement can be readily satisfied, provided that all the 
matrix elements of K are even integers, though. 

By varying Eq.(l) with respect to the Lagrange multi- 
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pliers a° a , one obtains a set of constraints 
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which determine the average values of the effective fields 
b a = B— < V x a Q > experienced by the CDF a-species. 

In the FQHE states viewed as the CDF IQHE, each 
of the CDF species occupies an integer number v a = 
2-Kpajba = m a of the effective LLs. The total electronic 
filling factor is then given by the expression 
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where fa = diag[mi, . . . , rrijv]. 

Integrating the CDFs out in the standard manner, one 
obtains a quadratic Lagrangian for the vector fields 
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where Hfj(u>, q) is the CDF polarization operator. 

Next, by eliminating all the statistical fields, one de- 
rives the RPA-like formula for the physical electromag- 
netic response function x'^iq) = IT"- 1 (q) + 2irKe t j k q k /q 2 
where q k — (o;,q). Quantized values of the Hall conduc- 
tivity corresponding to the putative FQHE plateaus are 
given by the formula 
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where a H = (2tt / 'uj)ImJ\. xy \ u )q _,,o is a tensor of the CDF 
Hall conductivities. 

As one important example of this general construction, 
attaching two units of the a- type flux (J* < V x a Q >= 
±47r) to the CDFs of the same type is equivalent to choos- 
ing K = ±diag[2, . . . , 2], in which case Eq.(3) yields 
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Notably, the overall Hall conductivity (6) is given by a 
"parallel" combination of the conductivities of the indi- 
vidual species (each of which is, in turn, given by a "se- 
ries" combination of the responses to the physical electro- 
magnetic A and the corresponding statistical field a a ). 
This composition rule should be contrasted against the 
naive one, a xy = J2a *W( 2 J2a v a+ 1 ) (see, e.g., the last 
reference in Ref «£) , which would have resulted from a se- 
ries connection between the response to A and a parallel 



combination of all the statistical fields (or a single field 
that couples symmetrically to all the fermion species). 

A more general case of the diagonal matrix K = 
diag[2p\, . . . , 2pjv] gives rise to a formula similar to 
Eq. (6) (with the factor of 2 replaced by 2p a in the denom- 
inator of the a-term). Furthermore, any "entangled" way 
of attaching the fluxes described by a non-diagonal K- 
matrix yields an expression different from Eq.(6). Such 
alternative choices of the iif-matrix (which we do not 
consider in this work) would be physically appropriate if 
different CDF species formed mutually coherent states. 
In the spin-polarized (N = 2) case, one example of this 
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sort is provided by the matrix K = 
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which has 



been previously discussed in the context of the double- 
layer v — 1/2 system. 

The 5'J7(4)-symmetry of Eq.(l) gets lowered in the 
presence of various symmetry-breaking terms of the form 
SL = Yla 3 Ir ^a^a/3*/3 which C an be of both, single- 
particle and many-body, nature (observe that these terms 
retain their form after the statistical transformation, if 
the matrix A is diagonal) . 

At the mean-field (Hartree-Fock) level, the list of such 
terms includes the (exchange-enhanced) Zeeman term 
(Ax = Ezl ® 1 ® a z ) and a parity-odd mass term 
(Am = Afjj eg) 1 eg) 1 or Aa z Cg> 1 eg) a z ), where the first, 
second, and third factors refer to the sublattice, valley, 
and spin subspaces, correspondingly. These two terms 
lift, respectively, the spin and valley degeneracies of the 
th LL, thereby splitting it into four individual sub-levels. 
Notice that the parity-even mass terms (~ a z Cg> a z Cg> 1 
or a z Cg) a z Cg> a z ) which have been previously discussed in 
the context of spin-orbit coupling and Spin Hall Effect 
in graphene can only split the th LL into two sub-levels 
(with the Zeeman term present). 

The symmetry-breaking terms appear to be instrumen- 
tal for describing, e.g., the plateau transitions — > ±1, in 
which case all the four sub-levels of the th LL are fully 
spin- and valley-resolved, as suggested by the strong-field 
(B > 20T) data of Ref£. 

The nearby FQHE states at \v\ < 1 can then be con- 
structed with the use of single-component (N — 1) CDFs 
which occupy an integer number of the effective LLs. 
Naturally, these states fall into the standard Jain's se- 
ries converging towards v = ±1/2 
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where m — 1,2, .. . and the overall ± sign is not cor- 
related with that in the definition of v^. Similar frac- 
tions can occur near v = ±3/2, thereby giving rise to the 
FQHE plateaus at a xy — ±(1 + v^). Since the single- 
particle CDF states are non-degenerate, a relative stabil- 
ity of the even vs odd-numerator fractions (7) is not an 
issue (cf. with the discussion in the last of Ref. 9 ). 

In the case of a residual SU(2) degeneracy of either 
spin or valley origin, the number of relevant CDF species 
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becomes N = 2. Conceivably, such a situation can occur 
at the — ► ±2 plateau transitions (where, say, v\ R = 

1/2 or i/^ = 1/2, depending on the relative magnitude 
of E z and A). 

The data of Ref. 5 suggest that at moderately strong 
fields (10T < B < 20T) the spin degeneracy gets lifted 
first (at least, at the n = ±1 LLs). In this scenario, 
the residual valley degeneracy gives rise to a series of 
valley-unpolarized IQHE states of the N = 2 CDFs which 
converge towards v = ±1 and correspond to the plateaus 
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where is defined in Eq.(7). 

In Eq. (8) , we took into account the fact that the num- 
bers of occupied (spin-polarized) effective LLs for the L— 
and R— type CDFs differ by one as a result of the spec- 
tral anomaly at the 0"* CDF LL (the i?-type states re- 
side at the energy E — A, whereas the L-type ones are 
at E — — A). As a result, the partial Hall conductivities 
of the R- and L-species as functions of chemical poten- 
tial obey the relation cr^ y (—u) = — <t£,(/x), although their 
sum a^y + a^y is, of course, an odd function of fi. It is 
worth noting that, from a formal standpoint, the anoma- 
lous IQHE observed in Refsi 2 - has the very same origin. 

Lastly, SU (4)-invariant spin- and valley-unpolarized 
states would be described in terms of N = 4 CDFs which 
provide a mean- field picture of the — 2 — > 2 plateau tran- 
sition in terms of the half-filled th LL {v^ R = 1/2) 
which is appropriate at relatively weak fields (B < 10T), 
according to the data of Ref£. 

Incompressible spin- and valley-unpolarized N = 4 
CDF states would then correspond to the plateaus 
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Notably, the series (9) includes (pseudo)spin-singlet 
states at v = 2/3 and 2/5, thus providing a possible 
CDF picture of the exact ground states found at these 
filling factors in the recent numerical studies^. 

By analogy with the conventional 2DEG4, we conjec- 
ture that the parent CDF states at v^ N=1 ^ = k - 1/2, 
(k = -1,0,1,2), i/( JV = 2 ) = ±1 and = for 

N = 1,2, and 4, respectively, behave as compressible 
"CDF metals" characterized by the presence of a Fermi 
surface of radius k* F = (2VB/N) 1 / 2 . 

The mean-field CDF dispersion relation remains lin- 
ear and the effective CDF velocity determined by the 
strength of the Coulomb interaction, v* F ~ gvpN 1 / 2 , is 
comparable to vp for g and N of order one. Due to their 
inherited Dirac kinematics, the Subnikov-de-Haas oscil- 
lations of the CDF resistivity at small deviations from 
the compressible fractions v^ N ^ are expected to show the 
same Berry phase of 7r as that of the original Dirac quasi- 
particles in weak fields 2 . 

The CDF Fermi energy E* F = v* F k* F ~ gvpB 1 / 2 ap- 
pears to be of the same order as the distance E\ — Eq — 



vp{2B) x l 2 between the th and ±l th LLs, suggesting that 
in graphene the LL mixing effects are potentially more 
important than in the conventional 2DEG where they 
get suppressed with increasing field. Moreover, the LL 
mixing becomes stronger with an increasing number n of 
the occupied electronic LLs, as the distance between the 
adjacent levels decreases as ~ |n| -1 / 2 . It is, therefore, 
likely that most favorable for the formation of the parent 
CDF metals and their incompressible descendants is the 
th LL (cf. with the conclusions drawn in Refs^ where 
the LL mixing was neglected from the outset). 

In the CDF IQHE states (7,8,9), the energy gaps for 
well separated particle-hole excitations 
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scale as ~ y/B/m for large m, which dependence is 
similar to that found in the conventional case of "non- 
relativistic" composite fermions where the effective mass 
varies as ~ B 1 ! 2 (see Ref.=). In contrast, for small m the 
true lowest energy excitations are likely to be represented 
by pairs of spin/valley (anti)skyrmions^. 

Despite the general possibility for the compressible 
states to emerge at any of the aforementioned fractions 
relative stability of these states is strongly de- 
pendent on the number N of the CDF species involved. 
In order to proceed with the stability analysis one has 
to go beyond the mean-field picture by including fluctu- 
ations of the statistical fields a Q controlled by the CDF 
polarization operator. 

In the regime where typical CDF energies and mo- 
menta are small compared to k* F and E F , respectively, 
tlij(q) is similar to that of a " non-relativistic" system 
with the same kp and vp. In particular, its transverse 
(with respect to the transferred momentum q) compo- 
nent n- L (u;,q) = q x n x q/q 2 = aq 2 + ibuu/q, where 
a ~ vp/kp and b ~ k* F , accounts for the Landau dia- 
magnetism and damping in the CDF metal. 

For N > 1 the CDF interactions are dominated by 
N — 1 linear combinations of the transverse components 
of the statistical fields which are orthogonal to the " in- 
phase" mode a a- Unlike the latter, these combina- 
tions are not affected by the unscreened Coulomb inter- 
actions, and the effective coupling between different CDF 
species (here V q = g/q, (v±v' ± ) = (vv') - (vq)(v'q)/q 2 ) 
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is always attractive in the Cooper channel (v = — v'). 
For N = 2 this interaction can facilitate the onset of s- 
wave valley-singlet pairing 10 . Moreover, for N = 4 there 
exists a possibility of more exotic (spin-valley coupled) 
patterns of the SU (4)-symmetry breaking. 

By contrast, for JV = 1 the effective interaction is re- 
pulsive in the Cooper channel, as it is between any CDF 
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species of the same kind for iV > 1, 

Although there is still a possibility of p-wave pairing be- 
tween the like CDFs, this potential instability (which is 
also present for the v = 1/2 state in the conventional 
2DEG) tends to be much weaker—. 

Since the inherent pairing instabilities make the N > 
1 CDF metals prone to becoming incompressible paired 
states, it is conceivable that the compressible states at 
the filling factors v^ 1 ^ should, in general, be less robust 
than those at z/ 1 ). Likewise, the chances of observing 
the novel series (8) and (9) might be rather limited, as 
compared to the standard one given by Eq.(7). 

The CDF metals would also be highly sensitive to dis- 
order. In the presence of potential (short-range) impu- 
rities of density pi, the CDFs experience elastic scat- 
tering off of an effective random magnetic field whose 
vector potential is described by the Gaussian variance 
< Ai(q)A,-(-q) >= 167r 2 pt(<% - <Mj/q 2 )/q 

A transport rate for the CDF a-species can be esti- 
mated as r* ~ Eppi/\p a \, and the above analysis (in- 
cluding the role of the symmetry-breaking terms) per- 
tains to the regime where E z , A, vpB 1 ' 2 /\m\ > T*. By 
the same token, disorder makes it more difficult to resolve 
metallic states at fractions v ~ l/2p with p > 1. 

Evaluating the longitudinal conductivity of the CDF 
a-species as a xx « max[\p a \/ pi, 1], one obtains a rough 
estimate for the physical conductivity of the CDF metals 

1 N 

^ K 7HEn -.^ ( 13 ) 
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which dependence should be contrasted with that at zero 
field (in the case of Coulomb impurities, the latter is 

proportional to the total electron density X^Pq™)- In- 
terestingly enough, in the experiment of Ref^. the con- 
ductivity at the v — plateau was found to be of order 
a xx ~ 0.6, possibly suggesting a precursor of the forma- 
tion of the N = 4 CDF metal at weak fields. 

Also, by analogy with the situation in the conventional 
2DEG 12 , we predict that the conductivity of the CDF 
metals is going to be temperature dependent due to quan- 
tum interference corrections which dominate over weak- 
(anti)localization ones and behave as 

r 2 r 

Sa xx oc -In a* x In- or -In — (14) 



for N — 1 and N > 1, respectively, thus allowing one to 
discriminate (in principal) between the single- and multi- 
component CDF metals. 

Furthermore, despite the ostensibly Fermi-liquid-like 
properties of the CDF metals, the electron spectral func- 
tion 7mG(p, e) exhibits a distinctly non- Fermi- liquid be- 
havior. Repeating the calculations carried out in the case 
of the conventional 2DEG^, we find a tunneling I — V 
characteristics of the CDF metal 

I(V) cx exp[-Const(E* F /V) v ] (15) 

where n = 1 for N = 1 and 1/2 for N > 1. 

To conclude, in addition to a wealth of its other re- 
markable properties, graphene provides a natural venue 
for the merger between the notions of pseudo-relativistic 
quasiparticles and statistical transformation from elec- 
trons to CDFs. We predict that compressible CDF states 
are most likely to be observed at the \Au\ — 1 plateau 
transitions (e.g., — > 1) between the fully resolved sub- 
levels of the th LL, since a stronger LL mixing makes it 
more difficult for such states to form at the \n\ ^ LLs. 

In contrast, the would-be CDF metals associated with 
the |Az/| =2,4 transitions are generally more fragile due 
to their propensity towards pairing, which can drive these 
states incompressible. Moreover, we predict that the in- 
compressible CDF IQHE states can occur at both, the 
standard (7) as well as novel (8,9), fractions. 

As far as the practical possibility of testing these pre- 
dictions is concerned, the detrimental effect of disorder 
calls for performing experiments of Refs.— & in samples 
of substantially higher mobility. The anticipated experi- 
mental signatures of the CDF metals can then be probed 
by such well-established techniques as bulk tunneling, 
acoustic wave propagation, magnetic focusing, and other 
geometric resonances^. 

For one, if a compressible v = state were indeed 
to occur at B < 10T, its CDF excitations would be 
amenable to conventional electrostatic gating. This pre- 
diction should be contrasted with such a hallmark of the 
zero-field Dirac kinematics as the celebrated Klein's para- 
dox that would hinder any possibility of electrostatic con- 
finement of the electronic Dirac excitations with vanish- 
ing Fermi momentum 2 . 
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